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ABSTRACT 


Weiner and Keast observed that in an upward-refracting atmosphere, the relative sound pres- 
sure level versus range follows a characteristic “step” function. The observed step function has 
recently been predicted qualitatively and quantitatively by including the effects of small-scale 
turbulence in a parabolic equation (PE) calculation. [Gilbert et al ., J. Acoust. Soc. Am. 87 , 
2428-2437 (1990)]. The present paper compares the PE results to single-scattering calcula- 
tions based on the distorted-wave Born approximation (DWBA). The purpose is to obtain a 
better understanding of the physical mechanisms that produce the step-function. The PE cal- 
culations and DWBA calculations are compared to each other and to the data of Weiner and 
Keast for upwind propagation (strong upward refraction) and crosswind propagation (weak up- 
ward refraction) at frequencies of 424 Hz and 848 Hz. The DWBA calculations, which include 
only single scattering from turbulence, agree with the PE calculations and with the data in all 
cases except for upwind propagation at 848 Hz. Consequently, it appears that in all cases ex- 
cept one, the observed step function can be understood in terms of single scattering from an 
upward-refracted “skywave” into the refractive shadow zone. For upwind propagation at 848 
Hz, the DWBA calculation gives levels in the shadow zone that are much below both the PE 
and the data. 


INTRODUCTION 

Weiner and Keast 1 and others 2,3 have observed that for sound propagation in an upward- 
refracting atmosphere, the relative sound-pressure level 4 versus range can be represented as a 
“step function” (Fig. 1). Recently the observed step function has been predicted qualitatively 
and quantitatively by parabolic equation (PE) calculations that include the effects of small- 
scale turbulence. 5 

Figure 2 shows gray-scale plots of the PE calculation without turbulence and with tur- 
bulence. The upward-refracted wave is called the “skywave.” 6 In the plots with turbulence 
the skywave is still present although it has been noticeably modified by turbulence. For a 
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source and receiver near the ground, the average relative sound pressure level inside the sky- 
wave (region 1 in Fig. 1) is approximately 0 dB (spherical spreading) with or without turbu- 
lence. However, the region below the skywave (region 3 in Fig. 1) is dramatically affected by 
turbulence. Without turbulence a deep shadow zone is predicted by the PE calculation. With 
turbulence, sound is scattered from the skywave into the shadow zone, producing a relative 
sound pressure level that is fairly uniform on the average. The region between the skywave 
and shadow zone (region 2 in Fig. 1) is a transition region. The horizontal extent of region 
2 is a strong function of the strength of upward refraction. It is evident that, for a gray-scale 
plot with turbulence in Fig. 2, a horizontal “cut” through the plot at a particular receiver 
height will give a step function. 

Although the gray-scale plots of the PE calculations give a good qualitative picture for un- 
derstanding the step function, the PE calculations do not allow a simple physical explanation 
of the observed quantitative behavior of the relative sound pressure level versus range. For ex- 
ample, what causes the observed constant relative sound pressure level (spherical spreading) 
at long ranges? In the present paper a simpler calculation is presented which is based on sin- 
gle scattering out of the upward- refracted skywave. The simpler calculation, which uses the 
distorted- wave Born approximation (DWBA), 7 is compared to the PE calculation and to the 
data of Weiner and Keast. The objective is to gain insight into the physical mechanisms that 
produce the observed step function. 


I. THEORY 
A. Atmospheric model 

We want to compare the DWBA calculations to the PE calculations reported in Ref. 5. 
Hence we use the same atmospheric model as in Ref. 5 and assume that the effects of tur- 
bulence can be adequately represented by small fluctuations in the index of refraction. The 
total index of refraction is thus written as a steady deterministic part n d (R) plus a fluctuating 
stochastic part fi(R,t) where R — (x,y,z) and t is time. With this approximation for turbu- 
lence, the wavenumber is given by 


k(R,t) = k 0 [n d (R) + fi(R,t)] , (1) 

where k 0 is a reference wavenumber, n d = 1, and fi • Cl. The deterministic part of the index 
of refraction n d is assumed to vary only with the height above the ground z. It was computed 
from a logarithmic sound-speed profile of the form 


( Co + a ln(z/d), z > z 0 
l Co + a ln(z 0 /d) z < z 0 


( 2 ) 


where Co = 340 m/s, z 0 = 0.01 m, and d = 6 x 10 -3 m. The refraction parameter a is —.5 
m/s for crosswind propagation (weak upward refraction) and -2.0 m/s for upwind propagation 
(strong upward refraction). The deterministic parameters were chosen to fit the short-range 
data of Weiner and Keast. 1 
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In the calculations reported here, the fluctuation fi is assumed at any instant of time to be 
a function of R = (x,y, z). In Ref. 5 the turbulence model was two-dimensional so that sound 
propagated only in the x — z plane. As will be shown, for computing the average sound pres- 
sure level using a single scattering approximation, the two atmospheric models are equivalent. 

In Ref. 5 the stochastic wavenumber in Eq. (1) was used directly to calculate “snapshots” 
of the acoustic field. Here, we want to compute average levels so we need the autocorrelation 
function for \i. The autocorrelation function is defined by 


C(S) = {fi(R + S)fi(R)) , (3) 

where {) denotes an ensemble average over many realizations of //. (We assume an ensemble 
average and time average are equivalent.) For small-scale turbulence near the ground, C(5) 
can be approximated by a Gaussian autocorrelation function of the form 

C(S ) = ^2 e -(S|//|+52//|+5?//2) 5 (4) 

where fi 0 is the root-mean-square value of //, and l x , l y , and / 2 are the correlation lengths in the 
x, y and 2 directions, respectively. In numerical calculations isotropic turbulence was assumed 
(L = ly = lz = 0 - The input values for fi 0 and l (fi 0 = 1.42 x 10 -3 , and / = 1.1 m) were taken 
from measurements reported by Daigle. 8 


B. Ground impedance model 

The ground was modeled as a flat, locally reacting plane with an angle- independent com- 
plex impedance. Impedance values were obtained from the empirical formulas of Delaney and 
Bazley 9 using an effective flow resistivity of 300 rayls/cm. The resulting impedance values 
were 7.19 + *8.20 and 5.17 + *5.57 at 424 Hz and 848 Hz, respectively. 


C. Distorted-wave Born approximation (DWBA) calculations 

We consider a point source with angular frequency u in a turbulent atmosphere. At a par- 
ticular instant in time the solution for a point source (Green’s function) in a turbulent atmo- 
sphere satisfies 


V 2 G'(R, R') + k 2 0 {n d + y) 2 G{R , R') = -AwS(R - fr) , (5) 


where R! is the source location, and R is the receiver location. In the absence of turbulence 
(fi = 0) the Green’s function Go is given by 


V 2 G 0 (R, #) + kln 2 d G 0 (R,R!) = -4t tS(R - Rf) 


(6) 
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where Go is a refracted wave (i.e., a distorted wave) if n<* varies with height. In this section we 
shall consider both an undistorted spherical wave and a wave distorted by upward refraction. 

An integral equation for G which goes to Gq in the absence of turbulence can be written as 


G(R,R') = G 0 (R,R') + i- ! Go{R,R!')Sk\R'')G{k',R')d 3 R" , (7) 

47T J 


where Sk 2 = &o( n d + l 1 ) 2 ~ ^o n d — 2/ifc, o , since = 1 and /< < 1, 

Equation (7) allows us to write the total solution, G, as the solution in the absence of tur- 
bulence, Go, plus an integral which gives the turbulent contribution. However, since the un- 
known G appears in the integral, Eq.(7) is as difficult to solve exactly as is the original differ- 
ential equation. When Sk 2 is “small enough” the full solution G that appears in the integral 
can be approximated by G 0 . The approximation G = Go is generally known as the “Born 
approximation”. When Go is a refracted wave the approximation is often called the “distorted- 
wave Born approximation” or “DWBA.” 7 

Writing the turbulent contribution as SG , using G = Go, and Sk 2 = 2/j.k q we have 

SG=^J G 0 (R, R")n(R")Go{R", R')d 3 R" . (8) 

We want to calculate the time average of |G| 2 , which we assume is the same as an ensemble 
average. Denoting the average value of \G\ 2 as (|G| 2 ) and assuming a random phase between G 
and SG we have 

<l<?l 2 > = {| Go \ 2 ) + <|«?| 2 ) , (9) 


where 


d & I 2 ) = / G'oA R'")GlA", #)<„(#'>(«")} (10) 

X Go(R, R")Go(R", R')d 3 R , "d 3 R" . 

Now (fi(R'")n(R")) = C(S) where C is the autocorrelation function and S = R ,n — R” . Trans- 
forming from the variables (R", R'") to the variables (R" , S) gives 

X Go(R, R")Go{R", R')d 3 s d 3 R" . 
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For the sake of illustration we first consider an undistorted spherical wave in free space, i.e., 
we take Gq(R,R") = exp(ifc&|./2 — R"\)I\R — /2"|. For |/?'| = 0 (source at the origin) and 
\R"\ » |S|, the Green’s function Gq(R" + S,R' = 0) can be approximated as 


G„(R" + S, K = 0) = 


(12) 


where the k = koh, and n = (R" /\R"\). Similarly the Green’s function Go(R, R" + S) can be 
approximated as 


.ik 0 \R-R"\ 

Go(R,R" + S) * =; — e 

| R-R”\ 


(13) 


where k' = k 0 h ', and h 1 = (R — R")/\(R — i?")|. With these approximations for the free-space 
Green’s functions, we have 


(I | 2 ) = ^ [ -4 ,-U— s C0)<PS<PR" 

Vl 1 ' 4 tt 2 J I R" \ 2 \R- R" |2 v ' 


We now define a scattering function a(q) as 


a(q) = J e i? S C(S)d 3 S , 


where q — k' — k. Then Eq. (14) becomes 


(14) 


(15) 


(| 8G | 2 ) = P- [ -J— <r(g) ■ 1 „ d*R" 

' 4tt 2 J I R" p kh, \R- R" 12 


(16) 


In Eq. (16) we can identify l/|.ft"| 2 as the sound intensity I tnc incident on the scattering vol- 
ume and 1/| R — R" | 2 as the scattered intensity I scat that reaches the receiver. Hence we can 
write 


(ISGf) = ■^JlUR")a(q)I, M (^')<PR" . (17) 

Equation (17) has a useful physical interpretation (see Fig. 4). The average intensity of the 
sound reaching the receiver from a particular volume of space is proportional to the product of 
the incident intensity reaching the volume, the scattering strength of the volume, and the scat- 
tered intensity. The Appendix gives an analytic result for Eq. (17) for small-angle scattering. 
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In order to take upward refraction into account we use the following prescription: In Eg. 
(17) we replace the incident intensity 7, nc = 1/|7?"| 2 and the scattered intensity 7, nc = 1/|7? — 
R" | 2 with Ii nc = \Gpe{R") | 2 and I scat = \Gpe(R,R")\ 2 , respectively, where Gpe denotes the 
Green’s function without turbulence computed using the parabolic equation method described 
in Refs. 5 and 10. Writing the integrals in Cartesian coordinates we have 


(I 6G | 2 ) = J | Gpe(R'') 2 | G pe (R,R ") | 2 

x exp [ i(q x S x + q y S y + q z S z )] C(S X , S y , S z )dx" dy” dz" dS x dS y dS z 


(18) 


Since Gpe is azimuthally symmetric and I72"! >> |5|, we neglect the y-dependence in Gpe and 
integrate exp (iq y S y ) over y and obtain the 6-function result obtained in the Appendix. (See 
Eqs. (A2) - (A4).) Also we set q x = 0, as in the Appendix, and obtain 


(I 1 2 > = J H / *"(* - *") I Gpe(R") PI G pe (R, R") p (ig 

xe iq * s >C{S x ,Q,S z )dS x dS z dx"dz" . 

Using the general Gaussian autocorrelation function in Eq.(4) for the integrations over S x and 
S z , we have 


(I SG P) = I X " {X - x ") I a PE (x'\ z") p (20) 

x | Gpe(x, Z] x", z") I 2 e- q * l *t 4 dx"dz" 

In the parabolic equation method the quantity actually solved for is ^(r, 2 ), where Gpp(r, z) — 
[exp(ifcor) / y/r ]^(r, z), and r = yjx 2 + y 2 . Since the integral is now two-dimensional, we can 
set y to zero and let r = x. Then in terms of ^(x, z) we have 


(I SG | 2 ) 


y 0 k 0 l x l z f T / // , /x ,2 1 lTf / 12 

— Y x — J I 2 ) II - x i z -< z ) I 

xe-W A dx"dz" . 


( 21 ) 


In the numerical calculations we assumed the turbulence to be isotropic with a correlation 
length /. Hence we have finally 


(I SG | 2 ) 


^ J | x, ** p| *(*-*»;*,*") P 

xe-&' 2 l 4 dx"dz" , 


where rr ;/ goes from zero to x, and z n goes from zero to infinity. 


(22) 
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II. NUMERICAL CALCULATIONS 

A. Comparison of DWBA calculations with PE calculations and with experiment 

In Fig. 3 the DWBA calculations and the PE calculations are compared to each other 
and to the data of Weiner and Keast. The data is for octave bands of random noise between 
300-600 Hz and 600-1200 Hz, respectively. In both the DWBA and the PE calculations, the 
frequency was taken to be where fi and / 2 are the lowest and highest frequencies, re- 

spectively, in the octave bands considered. Section I gives the parameters for the atmospheric 
model and the ground impedance model used in the calculations. Note that the DWBA calcu- 
lations and the data are for the average relative sound pressure level while the PE calculations 
are a “snapshot” of the relative sound- pressure level and not the average level. However, the 
trend in a particular PE calculation is generally fairly close to the average level predicted by 
the corresponding DWBA calculation. 

The DWBA calculations, which include only single scattering from turbulence, give a good 
approximation to the average PE levels in all cases except for upwind propagation at 848 Hz. 
For upwind propagation at 848 Hz, the DWBA prediction deep in the shadow zone is much 
below both the PE and the data. 


B. Discussion of numerical results 

Some of the features of the curves in Fig. 3 can be understood in a straightforward 
way using the DWBA calculations. The deterministic (no turbulence) part of the Green’s 
function is G 0 and the stochastic part due to turbulence is SG. Near the source (regions 1 and 
2 in Fig. 1) we have | Go | 2 ^ (| SG | 2 ) while at long ranges (region 3 in Fig.l) we have just the 
reverse. Consequently, near the source, the shape of a given curve for relative sound pressure 
level versus range is governed by the deterministic sound-speed profile so the level is essentially 
what one would obtain from a calculation without turbulence. 

Since we have | Go | 2< C (| SG | 2 ) at long range, the relative sound pressure level is due al- 
most entirely to scattering from turbulence. In order to understand the long-range behavior of 
the curves in Fig. 3 we must make a more detailed analysis than was required at short range. 
As shown in the Appendix, the contribution to the relative sound pressure level from turbu- 
lence scattering in free space (with no refraction) diverges as the logarithm of the range. We 
expect similar behavior even with upward refraction over a finite impedance plane. Consider 
the situation in Fig. 4 where we have a scattering volume with an incident intensity 7, nc and a 
scattered intensity I aca t- The sound intensity incident on the scattering volume is proportional 
to 1/r 2 where r is the horizontal range to the receiver. The scattering volume itself is propor- 
tional to r 3 . The scattered intensity reaching the receiver from the scattering volume, like the 
incident intensity, is proportional to 1/r 2 . For a fixed scattering angle, the average scattered 
intensity from the volume is thus proportional to (1/r 2 ) x (r 3 ) x (1/r 2 ) = 1/r. Hence, as 
shown in Eq. (A9), we expect the relative sound pressure level to increase as the logarithm 
of the horizontal range. This behavior at long range is seen in the DWBA calculation for 
crosswind propagation (weak upward refraction) at 424 Hz. When there is significant upward 
refraction the height of the scattering volume is not proportional to the range but increases 
more rapidly than linearly with range. As a result, the scattering angle is not fixed but 
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increases with increasing range. Since the scattered intensity is reduced at larger scattering 
angles, the relative sound pressure level versus range is “flattened” so that a nearly constant 
relative sound pressure level is reached at long ranges. Because of the flattening effect caused 
by an increasing scattering angle, a nearly constant relative sound pressure level is seen in the 
DWBA calculation for upwind propagation (strong upward refraction) at 424 Hz. A similar 
behavior is seen for crosswind propagation (weak upward refraction) at 848 Hz. The flattening 
effect with weak upward refraction at 848 Hz is apparently due to the greater sensitivity to the 
scattering angle at the higher frequency. 

The DWBA calculation for upwind propagation (strong upward refraction) at 848 Hz falls 
off in the shadow zone much more rapidly than the PE calculation and the data. The major 
computational difference between the two calculations is that the PE calculation includes mul- 
tiple scattering while the DWBA calculation does not. Hence the disagreement indicates that 
for upwind propagation at 848 Hz, multiple scattering is important. This interpretation is sup- 
ported by a gray-scale plot for this case which shows the skywave greatly modified by turbu- 
lence so that the approximation G = Go is not valid. 


III. SUMMARY AND CONCLUSIONS 

We have compared distorted-wave Born approximation (DWBA) calculations to parabolic 
equation (PE) calculations and to the data of Wiener and Keast. In all cases except one, the 
DWBA calculations, which include only single scattering, predicted the step-function behavior 
of the relative sound pressure level versus range seen in both the data and the PE calculations. 
The important conclusion to be reached is that, in the presence of upward refraction, single 
scattering can give a relative sound pressure level that does not diverge as the logarithm of the 
range but rather is nearly constant at long range. Hence, in all cases except one, the observed 
step function can be understood in terms of single scattering from an upward- refracted sky- 
wave. 

For upwind propagation (strong upward refraction) at 848 Hz, the DWBA calculation 
grossly underestimated both the data and the PE calculation. In this case, the single scat- 
tering approximation G = Go was not valid in the skywave. To accurately predict multiple 
scattering of sound into the shadow zone, one must have a good predictive model for sound 
propagation in the skywave itself. Hence, it would be valuable to have measurements not only 
for the sound scattered into the refractive shadow, but also for the sound field in the skywave. 
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APPENDIX: ANALYTIC EXPRESSION FOR FREE-SPACE SCATTERING FROM 

ANISOTROPIC TURBULENCE 

To see the general behavior of Eq. (17) we can consider weak small- angle scattering in free 
space . For weak small-angle scattering we can use the Born approximation and obtain an an- 
alytic result for anisotropic turbulence. 

For small angle scattering we can let 1/|R"| 2 = 1/x" 2 and 1/|R — R " | 2 = 1 /(x — x") 2 . We 
could integrate Eq. (17) directly using a particular autocorrelation function such as a Gaus- 
sian. However, to obtain a more general result that does not assume any particular autocorre- 
lation function, it is convenient to return to the form in Eq. (14) which is written in terms of 
the autocorrelation function C(S X , S y , S z ), 


<l iG l 2 > -£/ 


1 

2 


C(S X , Sy , S z ) 


p i(<lxS x -\-qySy+qzSz) 

(x - x") 2 

dx" dy" dz" dS x dSy dSz . 


For small angles we can approximate q x as 


(Al) 


1- = >* (i + -TTZ-y 


X — x" 


=k o 


[x"(x — x") 


Similarly, 


q y = k 0 


x 


[x"(x — X W )J 


We now consider the integral over z": 


(A2) 


(A3) 


[ + °° e iSzkolx/x"(x-x")]z" dz » _ 2ir S 
J — oo 


ko- 


x H (x - x") 
2n x"{x — x") 
ko x 


S(S Z ) . 


(A4) 


The integral over y" gives a similar result with S z replaced by S y . Inserting the results 
from integrating over z" and y" into Eq. (Al) and integrating over S z and we have 

(| SG | 2 ) = § / e’ , * 5l C’(S' x ,0, 0) dx" dS x . (A5) 

We could integrate over S x and define a special scattering function 
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(Too = f C(S X , 0, 0) exp (iq x S x )dS x . However, we are considering small-angle propagation. Hence 
q is almost perpendicular to the propagation direction and we can therefore set q x to zero. 
Thus, integrating in the region between the source and receiver we have 

(I SG | 2 ) = 


where 

Coo 

For anisotropic turbulence having a Gaussian autocorrelation function (See Eq.(4)) we ob- 
tain 


U 2 too tx 

/ C(S X , O, O) dS x / dx" 
X Z J-oo Jo 


k 2 

K o n 
— » 

x 


(A6) 


[°° C{S X , O, 0)dS x . 

J — oo 


(A 7) 


(| SG | 2 > = ^ 4 kl l./x . (A8) 

Thus for weak small-angle scattering in free space, the scattering due to turbulence falls off 
inversely with the range and depends on the correlation length only in the direction of propa- 
gation. Note that, written in terms of the relative sound pressure level (RSPL), the contribu- 
tion from turbulence scattering diverges as the logarithm of the range. 

RSPL = 101og w (x 2 (l SGI 2 )) (A9) 

= iOlogjo (V^o lx) + 101og 10 (x) . 
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RELATIVE SOUND 
PRESSURE LEVEL (dB) 



Fig. 1. Characteristic “step” function for the relative sound pressure level versus range for 
sound propagation in an upward-refracting atmosphere. (From Gilbert et al. 5 ) 
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Fig. 2. Gray-scale plots of relative sound pressure level versus height and horizontal range 
for a non-turbulent atmosphere and a turbulent atmosphere. The frequency is 424 Hz, 
and the source height is 3.7 m (12 ft). Parameters for the atmospheric model and ground 
impedance model are given in the text. (From Gilbert et al. 5 ) 
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Fig. 3. Relative sound pressure level versus horizontal range for a refracting, turbulent at- 
mosphere. The connected circles are the data of Weiner and Keast. 1 The solid lines are 
parabolic equation results from Ref. 5. The connected x’s are distorted-wave Born ap- 
proximation calculations. Parameters for the atmospheric model and ground impedance 
model are given in the text. 
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Fig. 4. Schematic representation of scattering from turbulence. The quantity /, nc is the 
average intensity incident on a particular scattering volume, and I scat is the average 
scattered intensity. The total scattered intensity is obtained by integrating over the vol- 
ume between the source and receiver. 
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